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We investigate a simple model of neutrino mass based on SU(3)c <8> SU(3)l ® U(l)x gauge unifi- 
cation. The Yukawa coupling of the model has automatic lepton-number symmetry which is broken 
only by the self-couplings of the Higgs boson. At tree level neutrino spectrum contains three Dirac 
fermions, one massless and two degenerate in mass. At the two loop-level, neutrinos obtain Majo- 
rana masses and correct the tree-level result which naturally gives rise to an inverted hierarchy mass 
pattern and interesting mixing which can fit the current data with minor fine-tuning. In another 
scenarios, one can pick the scales such that the loop-induced Majorana mass matrix is bigger than 
the Dirac one and thus reproduces the usual seesaw mechanism. 
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I. INTRODUCTION 

The recent experimental results of SuperKamiokande Collaboration 0, KamLAND and SNO Q confirm that 
neutrinos have tiny masses and oscillates. This implies that the Standard Model (SM) of SU(2) l ® U(l)y theory must 
be extended. 

The solar and atmospheric neutrino oscillations are now firmly established J] . The Am 2 values and mixing angles 
are known with fair accuracy 0, Q 

Am 2 tm = 2.4(l±°Ji) x lfT 3 eV 2 , 
Am 2 oi = 7.92(1 ± 0.09) x 10" 5 eV 2 , 
sin 2 23 = 0.44(l±g;£), sin 2 # 12 = 0.314(1^), 

sin 2 6»i3 = 0.9±H x 10~ 2 . (1) 
The tritium experiments provide an upper bound on the absolute value of neutrino mass 

m, < 2.2 eV (2) 

A more strict bound 



rrii < 0.6 eV 



was found from the analysis of the latest cosmological data 8] . 

Since the data only provide the information about difference in m 2 , the neutrino mass pattern can either be almost 
degenerate, or hierarchical. Among the hierarchical possibilities, there are two types: normal hierarchical or inverted 
hierarchical. In the literature, most of the models explore normal hierarchical case. 

In this paper we will explore a model which naturally gives rise to three pseudo-Dirac neutrinos with inverted 
hierarchical mass pattern. 
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Among the possible extensions of the SM, a curious choice are the 3-3-1 models which are based on the simplest 
non-Abelian extension of the SM group, namely, the SU(3)c <8> SU(3)l <g> U(l)x 0[H1 The reason why these models 
are appealing has been exposed in many recent publications 

The model requires that the number of fermion 
families be a multiple of the quark color in order to cancel anomalies, which suggests an interesting connection between 
the number of flavors and the strong color group. If one further uses the condition of QCD asymptotic freedom, which 
is valid only if the number of families of quarks is to be less than five, it follows that N is equal to 3. In addition, 
the third quark generation has to be different from the first two, so this leads to the possible explanation of why top 
quark is uncharacteristically heavy. 

There are two main versions of the 3-3-1 models as far as lepton sector is concerned. In the minimal version, the 
charge conjugation of the right-handed charged lepton for each generation is combined with the usual SU (2) L doublet 
left-handed leptons components to form an 5/7(3) triplet (t^, Z, Z c )z, . No extra leptons are needed and there we shall 
call such models minimal 3-3-1 models. There is no right-handed neutrino in its minimal version. Another version 
adds a left-hand antineutrino to each usual SU{2)l doublet left-handed lepton to form a triplet, i. e., (u, I, v c )l [To|. 
These left-handed anti-neutrinos serve the role of the charge conjugation of the usual right-handed neutrinos which 
are required in the usual seesaw mechanism. We therefore call such models right-handed neutrino models (RHNM) . It 
is this type of model that we shall explore in this manuscript. Its main feature is that it requires only a more economic 
Higgs sector for breaking the gauge symmetry and generating the fermion mass. Among the new gauge bosons of this 
model, the non-self-conjugated neutral boson X° can have promising signature in accelerator experiments and it can 
also be the source of neutrino oscillations [T^] • 

The explanation of the smallness of the neutrino masses and the profile of their mixing as required by recent 
experiments have been a great puzzle in particle physics. In the past several years a great amount of papers have 
been devoted to its solution (on the neutrino mass in the minimal 3-3-1 model, see Refs. PJE1G1). 

The most popular mechanism is of course the seesaw model with a few very heavy right-handed, SU(2)l singlet, 
neutrinos. This type of model requires a new very high scale of 10 12 GeV or higher. An alternative mechanism 
for generating small neutrino masses, which may not requires such high scale, is to do it only as one or multiloop 
radiative corrections. In the framework of SU(2) £g> U(l) model, a famous example is the so-called Zee Model and its 
generalizations |l6lll7j . In the framework of the minimal 3-3-1 model, this mechanism has been considered in [l4j 
based on the Zee type mechanism i. e. by introducing a scalar singlet. 

In this paper we shall explore the alternative RHNM in its minimal form. It is shown, with minimal Higgs sector, 
that the Yukawa sector has automatic lepton-number conservation which is broken in the Higgs sector. At tree level 
the neutrino spectrum contains three Dirac fermions, one massless and two degenerate in mass. At the two-loop 
level, very much like one of the Zee Models, with the help of lepton-number violating Higgs couplings, neutrinos 
obtain Majorana masses and correct the tree-level result. Since the Majorana masses involve a new physics (SU(3)l 
breaking) scale, depending on the size of the scale there are two scenarios possible. In the first one, the SU (3)l 
breaking scale is chosen to be very high and as a result the right-handed Majorana mass matrix is still very large, 
even though it is two loop-induced, compared with the Dirac mass. In this case, the usual seesaw mechanism still 
applies. In another scenario, the SU(3)l breaking scale is chosen to be not much higher than the weak scale; in 
that case the following interesting pattern arises. This naturally gives rise to an interesting inverted hierarchy mass 
pattern and interesting mixing which can fit the current data with some fine-tuning (to make the tree-level Dirac 
mass of order AmJ jm ). This radiative correction naturally occurs without introducing extra scalar singlet that was 
needed in the minimal model |14|. This scenario gives rise to a pseudo-Dirac neutrino mass pattern. There are many 
discussions of pseudo-Dirac neutrino mass pattern in the literature 0, 0, |2(| , however our scenarios is different from 
all of them as we will discuss. Throughout the paper we shall try to keep each sector minimal and see what kind of 
neutrino pattern is produced in general in the context of RHNM. We shall not implement by hand any extra texture 
in order to generate special pattern that can fit data. 

This paper is organized as follows. In Sec. Uwe review the 3-3-1 model with right-handed neutrinos and introduce 
the Higgs content and the Yukawa couplings. The conserved charges C and B are introduced and lepton-number 
violating couplings in the scalar sector are discussed and the general mass matrix is presented in Sec. 11111 while in 
Sec. II VI we derive the mass matrix through two- loop corrections. The main neutrino phenomena are are presented in 
Sec. [Vj Finally, the last section is devoted to our conclusions. 

II. THE 3-3-1 MODEL WITH RIGHT-HANDED NEUTRINOS 

To frame the context, it is appropriate to recall briefly some relevant features of the 3-3-1 model with right-handed 
neutrinos jlO|. In this model, the leptons are in triplets, in which the third member is a right-handed neutrino: 



SI = (ylM,Nl) T ~ (1, 3, -1/3), 1% ~ (1, 1, -1), 



(3) 
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where a = 1, 2, 3 is a family index. Here the right-handed neutrino is denoted by Nl = {vr) c ' ■ Note the fact that there 
are three generations of leptons is a peculiar consequence of anomaly cancellation as discussed in the introduction. 
This is an interesting plus to this type of models. The first two generations of quarks are in antitriplets while the 
third one is in a triplet and each charged left-handed fermion field has its right-handed counterpart transforming as 
a singlet of the SU(3) l group 

QaL = (d aL , -u aL ,D aL , ) T ~ (3, 3, 0), a = 1, 2, (4) 



Q 3 l = (u 3 L,d 3L , T L f ~ (3, 3, 1/3), T R ~ (3, 1, 2/3), (5) 
D a R ~ (3, 1, -1/3), u aR ~ (3,1,2/3), 

d aR ~ (3,1,-1/3), a =1,2, 3. (6) 

Note that the five quarks d a R and D a R have the same quantum number and so are the four quarks u a R and Tr. Their 
identity are defined only by the convention of the Yukawa couplings that we adopt as will be explained later. Also 
note that the third generation has different gauge content compared with the first two generations which is required 
by the anomaly cancellation. The electric charge operator is given in the form 

where X is the U(l) gauge charge, Aj are the SU(3)l gauge charge. The non-self-conjugated gauge bosons are defined 
as 

V2 W+ = Wl - iWl V2 Y~ = Wl - iWl 

V2X° = W^-iWl (8) 

where W l are the gauge boson associated A,. The physical neutral, self-conjugated gauge bosons associated with 
generator A3 and As and X, besides the photon, are again related to Z, Z' through the mixing angle <j). 

The gauge symmetry breaking and fermion mass generation can be achieved with just three SU(3)l triplets 

P = (Pt, Pl Ptf - (1,3,2/3), 

V = (Vi, fh, %°) T ~ (1,3,-1/3), (9) 
X = (Xi, X2". X") T - (1,3,-1/3). 

Note that the scalars r\ and \ have the same quantum numbers. By convention, wc define \ 1° b e the one with 
nonzero (X3) and breaks SU(3)l, therefore (77°) = by convention. The necessary VEVs are 

(p) = (0,u/V2, 0) T , (77) - (v/V2, 0, 0) T , 

( X ) - (0, 0,^/V2) T (10) 

In general, (x?) can also be nonzero, however, its effect is small and we shall ignore it in the following. Note that 
the identity of x is defined by convention to be the one will be responsible for SU(3)l breaking while p and r\ are 
responsible for SU(2)l breaking. The reason why two triplets are needed for SU(2)l breaking is because the three 
generations have different gauge charge, and one triplet is not enough to give fermion mass to all three generations 
(see below). 

The most general Yukawa Lagrangian as follows: 

L\ = KQzlTrx + h 2a pQ a LD l 3RX* + h.c. 

= ^i("3LXi + hhX2 + T L xl)T R + h 2a p(d aL Xi* - u aL xt + D aL xf)D f3R + h.c. 
£y = h ZaQsLUaRj] + h 4aa Q aL d aR ri* + h.c. 

= h 3a (u 3L ?]1 + d 3L r]- + f L r)l)u aR + h iaa (d aL rjl* - u aL i]+ + D aL T]^*)d aR + h.c. 
C P Y = h 5a Q 3L d aRP + h 6aa Q aL u aRP * + G ab fll b R p + F ab e 2jk (f L r(fE) bl (p*) k + h.c. 

= h 5a (u 3L pl + d 3L p° 2 + f L pl)d aR + h 6aa (d aL pi - u aL P2* + D aL p 3 )u aR 
+G ab [Vlp+ + l a L pl + N a L pt\l b R + 

+Fa b {mi C L ) b P^ - {N2) b pi]+l a L[{N2) b Pi {v C L ) b Pz] 

+N a L [{v C L ) h pl-{l C L ) h p- l ]} + ^c.. (11) 
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Note that, by convention, Tr is defined to be the one that couples to QzlX among the four quarks with the same 
quantum numbers, similarly, D^r, {(3 = 1,2), are defined to be the two quarks that couple to Q a LX among the five 
quark with the same quantum numbers. Note that G a b gives rise to charged lepton Dirac masses, while F a b, which is 
antisymmetric, gives rise to the Dirac masses for neutral leptons. 

The leptons have the Yukawa couplings only with the p Higgs boson. One can find a naive lepton number Ljv is 
violated only through the F a b coefficients while the rest of the whole Lagrangian, including the Yukawa couplings 
G a b, is the lepton-number conserving. Only the leptons carry the Ljv charge: Ljv(/^) = l,ijv(/£) = 1. Since 
phenomenologically, one requires F ao to be much smaller than G a b, it can be done in our context only by fine-tuning. 
The Ln allows us to claim that this fine-tuning is technically natural (in t'Hooft sense). We call Lm "naive" because it 
defines the lepton number z/^ = to be different from v^. Later, we will introduce another lepton number, L, with 
= L(vr) like the conventional lepton number, which will play an important role in our discussion on neutrino 
masses. The lepton Yukawa couplings needed in this work are presented in Fig. ^ 





FIG. 1: The necessary Yukawa couplings 



The VEV (x) breaks SU(3)l x U(1)x down to SU(2) L x U(l)y and gives masses of the exotic quarks as well as 
non-SM gauge bosons X,Y and Z'. The VEV (rj) gives mass for U3,d a quarks, while (p) gives mass for u a ,d^ and 
all ordinary leptons. The SM gauge bosons gain mass both from VEVs of 77 and p. 

After symmetry breaking the gauge bosons gain masses 

m 2 w = \g 2 (u 2 + v% M|. = i ff V+- 2 ), 

M\ = JsV+<A (12) 

W 4 and W 5 accidentally have the same mass. Ea. p2|) implies = u 2 +v 2 — 246 2 GeV 2 . In order to be consistent 
with the low energy phenomenology we have to assume that (x) 3> (p), (tj), such that mw Mx,My- The 
symmetry-breaking hierarchy gives us splitting on the bilepton masses | M x — My \ < . Since mw "C Mx , My , 
we can take Mx ~ My. The "wrong" muon decay limit 

r = y : - e \ v -^\ - feV < 1.2%, 95% cl 



T(p —> e UeUft) \M- 

gives My- > 230 GeV. From consideration of muon decay parameters, one has got the mass bound of the singly- 
charged bilepton of 440 GeV [2l|. With this mass scale, (x) - 800 GeV. 
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TABLE I: Nonzero lepton number L of fields in the 3-3-1 model with RH neutrinos. 



Fields 


N L l L 




Pi 


V$ 


Xi X 2 


D a L D/3L Tl Tr 


L 


-1 1 


1 


-2 


-2 


2 2 


2 2-2-2 



TABLE II: B and C charges for multiplets in the 3-3-1 model with RH neutrinos. 



Multiplet 


X 


n 


P 


Q3L 






d a R 


Tr D a R 


/ai 


laR 


B charge 











1 

-V- 


1 

-3- 


1 

3 


1 

3 


1 1 

3 3 








C charge 


4 
3 


'2 

3 


'2 

3 


3 


3 








-2 2 


1 

3 


1 



III. RADIATIVE CORRECTIONS TO NEUTRINO MASS 

The Yukawa couplings of Ea.ljlip possess extra global symmetries which are not broken by VEVs u,v,uj. From 
the Yukawa couplings, one can find the following lepton symmetry L as in Table [I] (only the fields with nonzero L is 
listed, all other fields have vanishing L). 

It is interesting that the exotic quarks also carry the lepton number. However, this L obviously does not commute 
with gauge symmetry. One can construct a new conserved charge C through L by making the linear combination 
L = x\z + y\$ + zX + CI where A3 and As are SU(3)l generators. One finds the following solution (see also p^ ) 

x = °' y = TI' 2 = °' and 

L = ^=\ 8 + £I (13) 
v3 

as in Table [n] Another useful conserved charge B is usual baryon number B = BI. 

Note that even though 77 and \ triplets have the same quantum number, they are distinguished already by our 
convention of Yukawa couplings and VEV's and, as a result, their lepton number assignments are quite different. rf( 
and X3 do not have lepton number L = 0, while 773 and Xi are bilepton L = 2. 

The lepton number C is, however, broken in the Higgs potential in general. The most general potential can then 
be written as the sum of the C conserving Vlnc an d £ violating Vlnv ( see also 

V(r), p, x) = V LNC (v, P, x) + V LNV (v, P, X), (I 4 ) 

where Vlnc{i1, P, X) is 

v L nc(v,p,x) = piv + v + pIp + p + pIx + x + >^i{v + v) 2 + M{p + p) 2 

+Hx + xf + (v + v)[^(p + p) + A 5 (x + x)] + \ 6 (p + p)(x + x) + x 7 (p + v)(v + p) 

+MX + V){V + X) + HP + X)(X + P) + {^e ijk niPjXh + h.c], (15) 

and 

V L nv{v,P,x) = jM 2 (x + V + V + x) + (V + X) [An {p + p) + A12 (r? + r?) + Ai 3 (x + x)] + h.c 

+Aio(x + ?? + V + x) 2 + [\uri + PP + X + h.c], (16) 

where overbars have been used to denote lepton-number violating couplings. The Higgs boson couplings necessary in 
this work are depicted in Fig. [21 

At tree level the neutrinos get Dirac masses from the Yukawa coupling F^Eij}. (/l) m {fi 1 ) b: '(P*) k - One can 
always assume that G a b is diagonal by convention and can always pick fermion phases so that the three coefficients 
F12, F13, F23 are all real. The resulting Dirac mass matrix is traceless and antisymmetric, and therefore has the mass 
pattern 0, — m„, m v . This is clearly not realistic. However, this pattern is severely changed by the quantum effect. In 
the base of (z/ e , v^Vt, N e , N^, N t )l, the most general mass matrix can be written as 
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Xk / \ Vi 



Pi 



FIG. 2: The necessary Higgs boson couplings 



M vN = 



f 









F\2 




M U /{ P2 ) 















— F13 


— -F23 





—Fyi 


— -Pi 3 






F\2 





— -F23 




Mn/(p2 


\F 13 


•P23 










Fl 3 \ 

F23 




/ 



(Pa 



(17) 



where My and M/v can arise from quantum correction. In particular, M v can be due to the loop-induced operator 

1 



Ou{M v ) ~ ^(/i/;,0(»7 + )V) J '> 
which is lepton number violating interaction, while Af/v is due to 

The Dirac masses can also receive quantum correction from the lepton number conserving operator 



O d (M d ) ~ ^{fifjK\ 



(18) 



(19) 



(20) 



In Ref. [24[, the effective dimension-five operators (O v , On) and {O4) were used to obtain the neutrino mass matrix. 
Choosing the free parameters in above operators (/, h and g) and taking v v = 10 2 GeV, v x = 10 3 GeV and A = 10 14 
GeV, one have got neutrino masses (mi, m 2 , ... me) in the range 10~ 5 -j- 1.7 eV. The authors argued that this set of 
parameters accommodates the solar and atmospheric oscillation data along with the LSND experiment altogether. 

It is well known that the neutrinos can get a mass through radiative mechanism [lfj . In the current model, neutrinos 
can get mass through two-loop radiative corrections, which is represented by the Feynman diagrams depicted in Fig|3] 

Since (77) <C (x) it is obvious that M v is small and negligible. The quantum corrections to the Dirac mass terms 
are also clearly smaller and negligible for the same reason [oc (rj) 2 , see a notice after Eg. 1)34(1] . However, radiative 
contributions to Mn can be very large and play a major role in determining the neutrino mass pattern in this model. 
The size of depends on the scale of SU(3)l breaking, (x) = and the dominant scale in the loop M in 

above equations. If uj/M is tuned to be very large, such as 10 12 GeV, then the M/v can be much larger than the tree 
level Dirac mass matrix and the usual seesaw mechanism will still apply. (In fact, in this case the loop corrections to 
the Dirac mass may also have to be taken into account). This scenario is more standard and requires very large ui/M 
which is not natural in our context. In the following we will concentrate more on the second, more natural, scenario 
in which ui/M is not so large. In that case, M/v can be considered a small correction to the dominant tree level Dirac 
mass. A curiously interesting neutrino mass pattern emerges. 



To finish this section, we mention that in Ref. |24j , the mass matrix obtained from dimension-five effective operators 
in 117(1 can give the possibility of explaining the LSND data as well as solar and atmospheric neutrino oscillation. On 
the other hand, depending on the temperature at which light sterile neutrinos thermalize, they can play an important 
role in big-bang nucleosynthesis (BBN) [25| or other aspects of cosmological problems 26] . All these problems should 
be further studied but it is out of the scope of the present work. 
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IV. TWO-LOOP CORRECTIONS 



Radiative correction stars only from the two- loop level. The most important two-loop contributions (to M/v) turn 
out to happen at the two-loop level. They are shown in Fig[3] 





X 

i X 




X >< 




p i 




p 1 • 


P :' 


i \p 




i \p 




f Fb£f c *) b F ab /«« 
a 


(/' !'" /•;., ./■ 


F»if c: y> F ab /«« 

b 




i X 


A 1X 


x i x 


P : 


»5 P t 

i \P 


P / 


' P i " ' * \ 



c d 



— ^Aii^ A 14 

p f 



if c yiG%ii a r G b u c ") h F ab /< 













^Aii 5 A 14 








p 


P-5 
P 





U C Y j F bd f G bci I- G?S e f 



FIG. 3: Two-loop contribution to neutrino mass matrix 

To calculate these diagrams, take [3] as an example, using vertices in Fig^ and FigEl the contribution of diagram 
is given by 



A d4k - p f P 'H + m fc) / i\ir*Mmp W + mJ jip f P 

{Ifabtiki' rii ^ § — \ l ) t bc e ^ L — cdtljt^R 



(2tt) 4 (2tt) 4 



= -2\ 14 Sf6^ 5 F ab F b * c F cd P R 
1 1 



fc 2 — TO 2 . 



(q 2 — TO 2 ) (fc 2 — to 2 ) (fc 2 — to 2 ) [(fc — q) 2 — to 
d 4 q d 4 k \ (-</) $ 



(2tt) 4 (2tt) 4 



(q 2 — TO 2 ) (fc 2 — TO 2 ) (fc 2 — TO 2 ) [(fc — q) 2 — TO 



q 2 — to? fc 2 — to 2 



Pr- 



(21) 



From l)2ip we see that only the neutrinos (i = j = 3), and not the leptons, received Majorana masses. Note that the 
integral in pip is finite. The contribution from both Figs. EJa) and b) (we call these F-type contributions) to the 
neutrino mass matrix is given 



M N{ai,d6) = -2^3,7i 5 Ai4^ 2 ^F ab F+F C(i A(6, C ) 

be 

= dij ^ FabFb c F cd A(b, c) 



be 



where it has been denoted a%j = —ISzi&zj^Xi/^o 2 and 

A(b,c) = [h{ml,m 2 c ) - h{m 2 c ,m 2 b )] = -A(c,b). 



(22) 



(23) 
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The factor 2 in a,y is due to summation over the SU(3) indexes. 
Here 



It is easily to check out that 



where 



ii(m 2 ,m 2 ) 



d A k d A q 1 

(2ir) 4 (2tt) 4 (fc 2 - ml) (fc 2 - m 2 ) 
4 l 



(q 2 - to 2 ) (q 2 - to 2 ) 



(24) 



[(fc — q) 2 — to 2 ] (fc 2 — m 2 ) 



FyiJ\ F13J2 

12 J\ F 2 3^ 

F13J2 F23J3 



M ( P ~ n 2 ( F 12 Ji F 23 J 3 I , (25) 



n 2 = 2^ 5 Ai 4 a; 2 , (26) 
J a = F 2 2 A(1, 2) + F 2 3 A(1, 3) - F| 3 A(2, 3), 



J 2 = F 2 2 A(l,2)+F 2 3 A(l,3)+Fi 3 A(2,3), 

J 3 = -^ 1 2 2 A(1,2)+F 1 2 3 A(1,3)+F 2 2 3 A(2,3). (27) 



We can approximate the integral by 



where M is the dominant mass scale in the loop: M ~ m p « m, ( . 

The contribution from FigsEJc) andEl^d) (we call these G-type contributions) to the neutrino mass matrix is given 

= a y -^[F Qb G bc (G+) + F rfh G bc (G+)] 

6,c 

d 4 k d 4 q fc 1 1 ^ 1 1 



X 



(2tt) 4 (2tt) 4 (fc 2 - to 2 ) (fc 2 - to 2 ) (fc 2 - m 2 ) (q 2 - to 2 ) (q 2 - to 2 ) [(Jfc - q) 2 - to 2 ] 

(29) 

Since TO b , to c <C m p ,m rp the loop integral depends on m b ,TO c only weakly. We see that the contribution is approxi- 
mately proportional to {FGG*) a d- Similarly, the contribution from FigsJS^e) and[3{f) to the neutrino mass matrix is 
given by 

MS = -a ij J2[ F a>>Gbc(Gt d ) + F db G bc (G+)] 

b,c 

d 4 k d A q \ k 1 1 4 1 



(2tt) 4 (2tt) 4 [ (fc 2 - TO 2 ) (fc 2 - TO 2 ) {q 2 - TO 2 ) {q 2 - TO 2 ) 

1 1 



[(fc — q) 2 — to 2 ] [(fc — q) 2 — to 2 ] 



(30) 



Note that the minus sign in Eq. 1(50) 1 is again due to summation over the SU(3) indexes. DiagramsEIc) -Elf) give a 
total contribution 
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= M^+M^=a t3 Y,[F ab G bc (Gt d ) + F db G bc (Gt a )] 

b,c 

d 4 k d i q f yfc ^ 1 1 



(2tt) 4 (2tt) 4 [ (fc 2 - m 2 ) (q 2 - to 2 ) {q 2 - to 2 ) (fc 2 - to 2 ) 
11 11 



(fc 2 — to 2 ) [(fc — q) 2 — to 2 ] [(fc — </) 2 — to 2 ] [(fc — q) 2 — to 2 
With the above mentioned approximation (that the integral being relatively insensitive to m b , m c ) we have 



(31) 



MjQ, * a 13 Y,[Fa b G bc {G + cd ) + F db G bc {G + ca )} (32) 

The contribution is again proportional to (FGG + ) T = G*G T F T = —G*G T F. This shows that two-loop contribu- 
tion as expected, is symmetric. Note that only coupling constant A14 contributes. In terms of the mass matrix, the 
G-type contribution gives 

A4 G) ~ pJ2 [Fa b G bc {G+ cd ) + F db G bc (G+ a )] , (33) 

b,c 



where 



1 \ 2 1 



P = n 2 [^) W (34) 



as before, n 2 = 2/i 5 Ai4W 2 . Two-loop contributions to M„ have the similar forms with just ay is replaced by n.^ = 
+25u5ij \14fA5v 2 and the mass of r\ in propagator is replaced by the mass of the \ boson. Note the plus sign on rib- 
and in this case is nonzero if i = j = 1. The VEV corresponds to the first component in the 77 triplet. 
Phcnomcnologically, it is necessary to fine-tune such that 

G ab » F ab , (35) 

because G ab is the charged lepton mass, while F ab is the neutrino Dirac mass. However, such fine-tuning is technically 
due to the protection Ln symmetry as discussed earlier. Therefore the F-type contributions in diagrams l^a) and 
EJb) are negligible and, hence Mm — Afjy . 

To look at M more closely , we can always choose a basis so that G ab ( and h 2a in the case of quarks) is diagonal. 
We have then 

(m^) = 2 P Y J F lb G bc G\ c = 2 P F ll G\ 1 = 0, (36) 

b,c 
b,c 

= 2pF 12 (\G 22 \ 2 -\G n \ 2 ) = (^4 G) ) 21 

Note that in the approximation (|35() . our result is similar to the one-loop radiative corrections ^?|]- Hence Eq.JSBJ) 
becomes 
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Noting that Gn ~ m e , G22 ~ fn^, G33 ~ m T , then the matrix in Eq. (|37|l has the form 




Denoting A = -^F 12 ,B = ^F% 3 , C = ^F 23 , then we can rewrite 

/ A B\ 

M d ~[ -A C (39) 

\ -a -c y 

Note that the relative size of relative to M d is controlled by the scale ratio (ui/M) 2 times some two-loop factor. As 
we states before, if the ratio (u>/M) 2 is chosen to be very large so as to overcome the two-loop suppression factor, the 
usual seesaw scenario can still apply. However, here we shall concentrate on the more interesting, and probably more 
natural, case when (oj/M) 2 times some two-loop factor is small such that Mjy can be considered a small perturbation 

to the Dirac mass M d . Next denoting r = ^ < 1, d = 2p F 12 (gft) = ^ P rAm 2 T , s = ^ P Bm 2 T , t = ^ P Cm 2 T 
then the neutrino mass matrix has the form 





1 









A 


B \ 






O 




—A 





c 


M vN = 








-B 


-C 








-A 


-B 





d 


s 




A 





-C 


d 





t 




\B 


C 





s 


t 


0/ 



(40) 



with d < s, t < A, B, C. We assume that u ~ u so u = ~ 175 GeV. 



V. PHENOMENOLOGY 



The interesting new physics compared with other 3-3-1 models is the neutrino physics. By our convention and 
from H39fl . it follows that Md is anti-Hermitian M~[ — —Aid, therefore its eigenvalues are imaginary and given 
0, ±iL, L = V 'A 2 + B 2 + C 2 . The eigenstates are 







( v e 








"3 J 







(41) 



where 



c jj(BC — iAL) -^(BC + iAL) 



'-(A 2 + C 2 ) jj{A 2 + C 2 ) = (f ,£ + ,f-) , (42) 



4 jt{AB + iCL) jj{AB — iCL) 



where L' = Ly/2(A 2 + C 2 ) and <pi are normalized eigenvectors of Md- The unitary matrix U3 diagonalizes Md, i.e. 

u 3 ul = ulu 3 = i, 

ulM d U 3 = D d = diag(0, -iL, iL) (43) 

We see that in the tree level we have three Dirac eigenstates. Two of them have degenerate eigenvalues L and the 
other one massless. It is easy to identify the mass splitting L as the value of measured atmospheric neutrino mass 
difference Am a t m . Therefore we require the parameter A, B, C to be of order Am a t m ~ 5 x 10~ 2 eV which is much 
smaller the charged lepton mass. This is of course part of the fine-tuning in fermion Yukawa couplings we need in this 
model. At loop level, this inverted spectrum is corrected by Mn, it will not only give rise to mass splitting between 
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the two degenerate Dirac states, it will also split each Dirac pairs into two non-degenerate Majorana states, resulting 
in the spectrum with six Majorana cigcnstates with four heavier ones and one light one and one remains massless. 
The existence of the massless Majorana state is a result of our approximation which gives M u n with zero determinant 
at the level of our approximation. We expect all the smaller (Majorana) mass splitting due to Mn should be of the 
order of Am = Am 2 ol / Am atm ~8x 10~ 4 eV. Here we are assuming that the solar oscillation is between the two 
heavier Majorana states. So for A, B, C of the same order of magnitude, 5 x 10~ 2 cV, we expect s, t to be of order 
8 x 1(T 4 eV. 

More specifically, with our loop correction d,s,t ^ 0, we can take them as perturbation and diagonalize the 6x6 
mass matrix. (Sec the appendix for more details). Note that, if we ignore CP violation, the mass matrix M u n is 
real, symmetric and therefore Hcrmitian. It can be diagonalized by an unitary matrix with real eigenvalues. When 
d = s = t = 0, the eigenvalues of M wN are L, 0, — L and eigenvectors of M vN are 

= (<f+, -i$+)/y/2, 

= i$-)/V2 (44) 



for eigenvalue L; 



*6 -¥-)/\/2 (45) 

for eigenvalue — L; and 

*T = (0, (46) 

for eigenvalue 0. In this basis, we can use Mm as perturbation and calculation the lowest order correction of d, s, t to 
the 6x6 mass matrix. The result (from the appendix) can be written as 




AM uN = A (47) 



where Aj arc the diagonal 2x2 blocks. Here we present only the diagonal blocks because they are the states with 
degenerate eigenvalues and therefore give the leading order corrections. The other off-diagonal blocks are nonzero, 
(they will be given in the appendix), but they only contribute at higher order. The Aj are 



A++/2 -A+_/2 
-A_+/2 A++/2 



(48) 



and 



where Ay are given by 



( o A°~ ) 



on 



A++ = —-(BCd — ACs + ABt), 

hi 



+ 



BC - iAL 

T 2 

AB + iCL 



L 2 



d + 
t, 



AC(A 2 + 2B 2 + C 2 ) + iB(C 2 - A 2 )L 
(A 2 + C 2 )L 2 



^oo 



-2A 



++ 



(50) 

(51) 
(52) 



They satisfy the properties A ++ = A , Ao- = A^ + , and Ay = A* { . The eigenvalues of matrix M„n are given 

by 



mi 



.2 = l + ^±l± 1 -VW7- 

A, 



(53) 



m 5 ,6 = -L + ^i±-v^A^F, 
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Note that A H is complex and 

IA+-I 2 - 



m 3 = 0, 

TO4 = — 2A ++ . 



{(A 2 + C 2 ) 2 [(B 2 C 2 + A 2 L 2 )d 2 + (A 2 B 2 + C 2 L 2 )t 2 



(54) 
(55) 



L 4 {A 2 + C 2 ) 2 

-2{A 2 + C 2 )ACdt] + (A 2 B 2 + C 2 L 2 )(B 2 C 2 + A 2 L 2 )s 2 

+ 2{A 2 + C 2 )\(B 2 C 2 + A 2 L 2 )ABds + (A 2 B 2 + C 2 L 2 )BCst]} . 



(56) 



Therefore y/\ A_| | 2 is nonanalytic in d, s, t. For example, in the simplified case of A = B = C = A4, we have 

L = MV3, A++ = l(d-s + t), |A+_| 2 = £(d 2 + s 2 + t 2 + ds + st - dt). 
In this case, the eigenmasses are given 

mi, 2 = A^\/3 + i(d - s + i) ± ^(d 2 + s 2 + t 2 + ds + st - dt), 
m 5 , 6 = -Al\/3 + - s + i) ± ^-\/(d 2 + s 2 + 1 2 + ds + st - dt), 

o 

to 3 = 0, 

2 

m 4 = --(d - s + t), 



(57) 
(58) 

(59) 
(60) 



The eigenvectors of to = and m = Aoo have a form 

1 

1 

For other eigenvectors with eigenvalues L, we get 

±( 1 

V2 V e ~ i4 

±( 1 

V2 V 



(* 3 + «*4) 
(* 3 -i* 4 ) 



_(* 1 + e- i ** 2 ) 



(61) 



V2 



(*i-e-*** 2 ) 



where is the phase of A + _ = | A H |e l< ^. The eigenvectors with eigenvalue — L are similar linear combination of ^5 

and * 6 . 

The mixing matrix in the basis v e ,y ll ,y T ,N\,N2, N 3 and v\, V2, ^3, V5, v§ is given by 



(62) 



/ Ve \ 










v\ 


V T 


= u 


vz 


N\ 


1/4 


N 2 






\nJ 




We/ 



where 



U = 



( u el u e2 

Un ± i Un i2 

Un 2 i Un 2 2 

\ Un 3 1 Ujy 32 



L 

_BL\/2 
L 

AL'V2 
L 

iCL'V2 
L 

■ BL's/2 



CL'^2. 

BL'V2 
L 

AL'V2 
L 

-iCL'^2 
L 



Url 
BL'V2 tt 



AL'^2 iAL'^2 tt tt 

— Z - U N 3 1 -Un 3 2 ) 



u e2 

Un 2 2 



(63) 
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TABLE III: Matrix elements 



11 , 




V e2 


BCk~ - {ALk+ 


11 1 


(A 2 + C 2 )k , 


11 „ 




ABk + + iCLk- 


u t2 


ABk_ + iCLn + 




-ALn+ - iBCn- 


Unx2 


-ALk- - iBCK + 


Un 2 i 


~i{A 2 + C 2 )K- 
~i{A 2 + C 2 )K+ 


Un 2 2 


Un 3 i 


CLk+ - iABk- 


Un 3 2 


CLk- - iABn + . 



where the matrix elements are given in Table HTT1 

Here we have denoted k± = 1 ± e~ 1 ^ . Note that C/jv x i = —iU e 2, Un 2 i — —iU^, Un 3 i = —iU T 2, Un x 2 = —iU e i, 
U N22 = -iU^i, U Ns2 = -iU Tl and ^ = 2VA 2 +C 2 . 

One can check out that this matrix is unitary 

UU* = U ] U = 1. (64) 
Here the inverted hierarchy neutrino mass spectrum is used and is shown in Fig0] 



<"2 

mi 






solar 


m 6 





m$ 


atmospheric 


m 4 





FIG. 4: The inverted hierarchy neutrino mass spectrum, showing the usual solar and atmospheric mass differences 



The survival probability is given by |27| (in the extremely relativistic limit) 

Pvi = |<^(o)h)WI 2 = El^^W^I 

( lux \ 

X COS tpil'aB , (65) 

where L a p = and tp u , a p is the phase of Ui a U*, a U*pUi' p , with Am^ =m% — nip. 

In the literature, usually the only simplest two component neutrino mixing result is used in the analysis. However, 
in our case of six light Majorana neutrinos, the analysis is much more complicated and has not been fully explored in 
the literature. Even the oscillation formula for three flavor which is available in the literature is far from useful here. 

For our special neutrino spectrum, however, we can make some approximation and get some nice result. Since the 
small iS.m 2 sol is mainly between v e and and the larger Am^ tm is mainly between v r and u^, it is reasonable to 
assume that v r is mainly contained in v% and v± which are very light, while v e and are mainly contained in v\, ^2, 
1/5, i>q which are almost degenerate in mass. In any case, this is almost the only reasonable guess for any theory with 
an inverted hierarchy neutrino mass spectrum. We shall assume this here. 

In this case we can get an approximate result for P^ T in the vacuum, by considering only the transition between 
light mass eigenstates ^3 and 1^4 and those of heavier v\ , V2 , ^5 , ■ Since the other transitions involving states whose 

1 
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mass splitting are too small for atmospheric range oscillation. Therefore 

P^t ~ E U ^U* a U* a U Ta + Y U Hf3 U T/3 U ll(3 U Tf3 
a=3,4 0=1,2,5,6 

+ E E \ u ^K a u; f3 u Tf3 \ 

a=3,4 0=1,2,5,6 



X COS ( - IffiTap) j (66) 

where U is given in Ea. (|63|l . for example, U^i = U ^ = — • >K+ , while C/ r 3 = ^Aj- £ is in the range of atmospheric 
neutrino oscillation p^. 

Substituting matrix elements into Eq. 1)66(1 we see that the first term in (|66[) is given by 

Y u, ia u; a u; a u Ta = ^L, (67) 

a=3,4 

and the second term in (|66|l has the form 

Y u^u; p u; p u TP = -^[a 2 b 2 + c 2 l 2 

0=1,2,5,6 

+ cos 2 cj)(A 2 B 2 ~ C 2 L 2 ) - ABCL sin 2<j>] (68) 

Let us calculate the third term 

Y V m U* ra = (69) 

a=3,4 



2. u m u rf> - — jofjf — - -jT- ( 7 °) 

0=1,2,5,6 v y 

Here two factors are all real. Therefore ^ TO (3 can be ignored. Hence 

2^ 2^ \U t , a U Ta U^Urp\ COS I If^ap J = — - COS I I . (71) 

a=3,4 0=1,2,5,6 \ a(3 J \ atm J 



We get then transition probability 



P^ T = -L [3A 2 B 2 + C 2 L 2 + cos 2 0(A 2 B 2 - C 2 L 2 ) - ABC 'sin 20 - 4,l 2 5 2 cos ( 



AL 4 ' " \L 

where 



(72) 



r _l _ ( TO 1,2,5.6 ^ m l,4) i 2 

Qtm " 4^p ~4^p lMj 



For is e , ^ T oscillation we have 



P eT « 53 U ea U* Ta Ul a U Ta + Y U e pU* Tp U* efs U Tfi + E I^Wa^^l COs(|p-¥>er«tf) (74) 

a=3,4 0=1,2, 5, 6 a=3,4 0=1,2,5,6 V Q / 

Using a similar approximation, we obtain 

V u u-iru - 32 ^ 2 (^ 2 + g2 ) 2 

a=3,4 ^ ' 

J] U e0 U^U* e0 U rf3 = -i— {8[i? 2 C 2 + A 2 L 2 + cos 2 4>{B 2 C 2 - A 2 L 2 ) 

0=1,2,5,6 * ' 

+ 2 ABCL sin <j>] [A 2 B 2 + C 2 L 2 + cos <j){A 2 B 2 - C 2 L 2 ) - 2 ABCL sin 0] 
+8[E 2 C 2 + A 2 i 2 - cos 2 0(B 2 C 2 - A 2 L 2 ) 

- 2 ABCL sin (/)}[A 2 B 2 + C 2 L 2 - cos <t)(A 2 B 2 - C 2 L 2 ) + 2ABCL sin </>]} , (75) 
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V u u* 

q=3,4 

E U:pUr, = - 2AC{A ^ +C2) . (76) 

(8=1,2,5,6 



These values again are real. Summing up, we get 



Per = jrj T -^^{2A 2 C 2 {A 2 +C 2 ) 2 + (A 2 B 2 +C 2 L 2 )(B 2 C 2 +A 2 L 2 



+ (A 2 B 2 - C 2 L 2 )[cos 2 <\>{B 2 C 2 - A 2 L 2 ) + ABCL sin 20] 
- ylBCL[(B 2 C 2 - A 2 L 2 ) sin 20 + 4ABCL sin 2 0]} - cos f |^ 

i \L a tn 

From (|77|l we see that to get sin 2 #13 w we just need C — 0, — or equivalently 

C = 0, d = s = 0. (78) 
In this limit, the atmospheric neutrino transition probability (|72() becomes 



^B 2 
L 4 

„ „:„2 



1 — COS 



27TX 



yl 2 S 2 



2- 



L 4 

(79) 



Eq. {7HJ) gets the usual form [2? 



by identification 



Thus 



P MT = J sin 2 2# atm sin 2 ( ) ( SI) ) 



A 2 £? 2 4A 2 £? 2 

sin 2 20 otm = 4-^- = ( ^ 2 + j32)2 (81) 



sin 2 20 Qtm = l => ,4 = 5. (82) 



To finish this step, we note that to get sin #13 ~ and sin 26 a tm — 1, one just need C — 0, A = B. 
It is much harder to make an approximate calculation for P eM since it involves 



T / 2irx \ 

Pep » 2J \ U e a U* a U*^U^\cOS f - <p ef i a f3 J , 

i 9 r. r \ a/3 / 



a/3=l,2,5,6 



(83) 

for a; is in the range of solar neutrino oscillation. It is hard to proceed analytically. However, we can make a numerical 
study of the possibility to make sin 2 9 so i ~ 0.3. Of course, the solar neutrino oscillation is mostly likely due to 
matter-induced MSW oscillation and our vacuum oscillation treatment is flawed, but it serve to illustrate that our 
model can easily fit the solar data also. A more detailed serious study of the six Majorana flavor (for this or other 
similar pseudo-Dirac models |18| ') will be needed. This and a study of the astrophysics constraint will be investigated 

in a future publication. 

Now we consider v e , transition and for the sake of shorthand denote q = \/|A_| | 2 . We have 



m 2 = I? ^ p- + \ +LA++ + Lq+ -A ++ <7, 



A 2 


<r 


+ r 


+ T + 


A 2 


4 + 


+ r 




A 2 


4- 


+ r 




A 2 

+¥ 


4- 







">i = ^ + ^— t- -V LA ++ - L <i- 2 A ++ 9 ' 

L 2 + + ~ - LA++ -Lq+ \&++q, 

ml = L 1 + ==± + ^ - LA++ + Lq - ±A ++ q (84) 
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Some manipulations give 

P ell {A,B,C, d,s,t) 



where 



4L 4 
1 



8L 4 



[{B 2 C 2 — A 2 L 2 ) cos (f> + 2ABCL sin 0] 
{cos <f) [{B 2 C 2 - A 2 L 2 ) cos + 2ABCL sin 0] 



L(A ++ +q) 



L(A ++ - q) 



sin cj> [(B 2 C 2 - A 2 L 2 ) sin - 2ABCL cos < 



X cos 



cos 



g(2L + A +H 



2m 2 
A++(2L-g) 
2m 2 



cos 



cos 



^ = 10- n [ey 2 l 



is the solar oscillation parameter. 

In the case C = 0, the v e Vy, transition probability becomes 



A ++ (2L + g) 
2m 2 
q(2L - A ++ 



2to 2 



(85) 



(86) 



A 2 

P efJ ,(A,B,0,d,s,t) = ^cos 

A 2 



L(A+ + + g) 



A ++ (2Z + g) 



L(A ++ -g) 



2m 2 



771 

A ++ (2L-q) 



2m 2 



cos 



g(2X + A 



++) 



2m? 



cos 



g(2L - A^ 

27T7 2 



Putting one more condition, d = (together with C = 0) we have L = y/ A 2 + B 2 and 



cos< 



S11K 



9 = 



At 


v / s 2 (A 2 +B 2 ) + 


A 2 t 2 ' 


sVA 2 + B 2 




^s 2 (A 2 +B 2 ) + 


A 2 t 2 ' 


By/s 2 {A 2 +B 2 ) 4 


- AH 2 



A 2 + B 2 



ABt 



A 2 + B 2 ' 

Substituting Eq.® and (JHH into JHZJl we get for s = 

P eiM (A,B,0,0,0,t) 



A 2 



A(A 2 + B 2 ) 8(A 2 + B 2 



1 + cos 



2 x 10 1 



ABt 



VA 2 + B 2 



(87) 



(88) 



(89) 



(90) 



In Fig. 0we plotted P efl (0.12, 0.1, 0, 0, 0, i[0.1eV]) for t run from to 10~ 8 or from to 10" 9 eV 

The figure shows that, if A and B are taken in order 0.1 eV (upper cosmological bound) we have solar neutrino 
data for t parameter. 

For the v e v^ transition in matter, ignoring magnetic field, we get typical terms similar to that in Ref. 19j. For 
general case, P a( g cannot be cast in compact form and may require numerical evaluation. Analysis of the matter effect 
is quite simple in the case of two generations and gives a result in good agreement with the current evaluation. 



17 



P_(emu) 

0.24 
0.22 



0.1B 
0.16 
0.14 

FIG. 5: P eM as function of t within A = B = 0.1 eV, C = d = s = 0. 



VI. SUMMARY AND CONCLUSIONS 

The basic motivation of this work is to study neutrino mass and mixing in the framework of the model based on 
SU(3)c <8> SU(3)l <g> XJ(l)x gauge group with right-handed neutrinos. 

The Higgs sector of this model contains the bilepton Higgs scalars with lepton number L = 2. Hence, the Yukawa 
coupling of the model has automatic lepton number symmetry which is broken only by the self-couplings of the Higgs 
boson. The interesting radiative mechanism for neutrino masses has been obtained. At the tree level, the neutrino 
spectrum contains three Dirac fermions, one massless and two degenerate in mass. At the two-loop level, neutrinos 
obtain Majorana masses and correct the tree-level result which naturally gives rise to the pseudo-Dirac mass differences 
and an inverted hierarchy mass pattern and interesting mixing which can fit the current data with minor fine-tuning. 

For the solar neutrino oscillations in matter, neglecting magnetic field, we have got the matter effects. Our analysis 
in the simpler case limited by two generations shows that the scheme gives appropriate consistency. 

In another scenarios, one can pick the scales such that the loop-induced Majorana mass matrix is bigger than the 
Dirac one and thus reproduce the usual seesaw mechanism. 

The complete analysis and study of the astrophysics constraint, the medium effects will be investigated in the future 
works. 
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APPENDIX A: DIAGONALIZATION OF THE MASS MATRIX 

Let us consider a free case 

Ho=[ m t J, (Al) 
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where M T = -M, M* = M and 



M 



A B 
-A C 
-B -C 



(A2) 



/O 0\ (0 d s\ 

The interaction is considered as perturbation V = I „ . I , with A = d t . Let us start with a simple case 

V U A / \a t J 



M 



A = then H = Hq + V = I ^ J ' ^ e searc ^ ei g envec t° rs m the form ^ e 
e = E^> e ■ Then we obtain an equation 



9 a 



satisfied the equation 



M 
-M 



<Pb 



E 



9a 

9>B 



Equivalently, 



M(p B = Eip A , 
-Mtp A = Eip B . 



(A3) 



From Eq. I|A3J| . we get —M 2 tpA = EMip B = E 2 lpa, or (— M 2 )ipA = E 2 ipA- This means that we need only to 
diagonalize M. Next, we consider a characteristic equation 



-A A B 
-A -A C i = (i 

—B -C -A 



(A4) 



or \{X 2 + A 2 + B 2 +C 2 ) = 0. Thus, we get three roots : A = 0, ±WA 2 + B 2 + C 2 . Let us denote L = VA 2 + B 2 +C 2 , 
then 



E z 



-X 2 = 0, -{±iLf = 0, L 2 , 



Thus we have three eigenvalues : +iL, 0, — iL. Let us choose if a to be M's eigenstate 



It is easily to get y 



-C 
B 

-A 
C 






A 


B 




-A 





C 


)(;)- 


-B 


-C 







chanj 


l& a 


sig 


;n of the eigenstate 



-c 



-A 



(A5) 



(A6) 



C 

— B | . Thus, we get a 



massless eigenstate | ) = j- \ —B . Now we look for other eigenstates 



.4 



A B 
-A C 

-B -C 



V = ±iL y 



or Ay + Bz = (±iL)x and — Ax + Cz = (±iLy) or ±iLx — Ay = Bz and Ax ± iLy = Cz. We get then, 

' BC — iAL A 2 + C 2 



AB + iCL ' AB + iCL ' 



1 



for A = iL, and 



-£?C - iAL 



A 2 + C 2 



-AB + iCL' -AB + iCL' 
for A = —iL. Thus the normalized eigenstates are given by 



(A7) 



(A8) 



(A9) 
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1 



y/2(A 2 + C 2 )(A 2 + B 2 + C 2 ) 
1 

V / 2(A 2 + C 2 )(A 2 +B 2 + d 2= ) 



{BC - iAL, A 2 + C 2 ,AB + iCL), 



(BC + iAL, A 2 + C 2 ,AB- iCL) 



(A10) 
(All) 



Now we have obtained three eigenstates | 4>o )■> I 4*+ )> 10-)- It is straightforward to check that they form 
an orthonormal basis. Note that, for an anti-Hcrmitian matrix like M, its eigenstates (with different E) are also 
orthogonal! 

Now, it is straightforward to construct the eigenstates of the true Hamiltonian H . For E ^ 0, Eq. l|A3Jl gives 



(f B = -—MtpA 

1. For E = L: ips = —-^M(fA- Taking tp A = (f) +1 then tps = — j^L(f> + = —i<f)+. Thus we have ^ 1 = ^= 



(A12) 



Similarly, let ipA = 4>-, then Lps = — j;(—iL)(f>- — i(f)-. Hence ^2 = ^75 



Note that (*i|* 2 ) 



.(1-1)=Q. 



2. For E = —L: ips = j^MipA- Analogously, we get \I> g = -7= 

3. For E = 0: a special case <5 3 = f ^° \ and <f 4 = f ^ 
Now we only need to rewrite the perturbation V in a new basis 

Vij = (*i|V|*j>, z,j = l,2,...,6, ^= (J a) 



and * 6 = ^ 



(A13) 



We also have A = £\ . \(j>i)Aij{<pj\. Sine are expressed through just three functions <p + , O , </>-, so we only need 
to work out the 3x3 matrix! Finally, the 6x6 matrix is 



/ 




-±A + _ 





73 A +° 


-^A ++ 


^A + - 


\ 






|A- 







-|A_ + 


-|A- 



































A 00 


7I A °+ 


-73 A °- 






-■§a ++ 


2-A+- 





-7f A +° 
73 A -° 


^A ++ 


-§A + _ 




V 




-|A- 







|A- 


/ 



To calculate the leading corrections, only need to diagonalize the 2x2 matrix. Note that Ay = Aij(d,s,t) are 





homogeneous function of order 1 for (d, s, t). The matrix ( ^ ^ j , gives AE — 0, Aoo + 0{g 2 ). which is analytic 
function. So, the energy spectra for E$, E4 are given 

+ O(g 2 ), E 4 = A 00 + 0(g 2 ) 



-A --A 

Next, the matrix ( 2 1 ^ + j 2 . + 

"2^ — y 2^ — 



E;i 



, has two eigenvalues given by 



(A14) 



AE= A ++ + A -± 



A++ - A__ 



IA, I 



(A15) 



which is, in general, nonanalytic!. For instance, f(d,s,t) = \/d 2 + s 2 + t 2 is not analytic at (d,s,t) = (0,0,0). Thus, 
we obtain 
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L 



A 



A__ 



£-5,6 — —L 



4 

A ++ + A_ 



± 



A^ 



A__ 



I A, I 



± 



A ++ - A_ 



0(g 2 ), 



0(g 2 ) 



(A16) 
(A17) 



The O(g) corrections are nonanalytic in general!. To get explicit results we calculate Ay. 
A++ = 0+A0+ = 2{A i+ c z )L z 

/ d s\ I BC — iAL 

x (BC + iAL, A 2 + C 2 , AB — iCL) I d t\\ A 2 + C 2 ] = ^(BCd- ACs + ABt). Noting 0_ = (f>* + , we have 

\ s £ / y AB + iCL 

A__=A ++ . 

Now we turn our attention to Aoo- Aoo 
ABt) = -2A ++ . 



d s 

6+A^o = ^(C,-B,A) \d0t 

s t 



A+_ = 0+A^_ = 



2( J 4 2 +C 2 )L ;i 





d 


















c 

-B 
A 

BC + iAL 
A 2 + C 2 
AB - iCL 



jjz(BCd — ACs 



= (^^)d 



AC(B 2 +L 2 )-iBL(C 2 -A 2 ) 
{A^+C 2 )^ 



We have A_+ = (A+_)*. Similarly A 



o+ 



0+A0 4 



A(L 2 -2B 2 )-iLBC 
^2(A 2 +C 2 )L 2 

A _ = #j-A0_ = 



and 

C(L 2 -1B 2 )-iABL 



y<2(A 2 +C 2 )L 2 

Noting that <fio is real, and <p + - 
calculate 



y / 2(A 2 +C 2 )L 2 



C(L 2 -2B 2 )+iABL 
^2(A 2 +C 2 )L 2 



A(L 2 -2B 2 )+iLBC 
^2(A 2 + C 2 )L 2 



y/2(A 2 +C 2 )L 2 



we have the following properties Aq_ = Aq + = A*_ = A + o- To complete, we 



IA+_| 



1 



{(A 2 + C 2 ) 2 [(B 2 C 2 + A 2 L 2 )d 2 + (A 2 B 2 + C 2 L 2 )t 2 



L 4 (A 2 + C 2 ) 2 

- 2(A 2 + C 2 )ACdt] + (A 2 B 2 + C 2 L 2 )(B 2 C 2 + A 2 L 2 )s 2 

+ 2{A 2 + C 2 )[{B 2 C 2 + A 2 L 2 )ABds + (A 2 B 2 + C 2 L 2 )BCst]} 



(A18) 
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